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Abstract. We present a new perturbative formulation of non-equilibrium thermal field theory, 
based upon non-homogeneous free propagators and time-dependent vertices. The resulting time- 
dependent diagrammatic perturbation series are free of pinch singularities without the need for 
quasi-particle approximation or effective resummation of finite widths. After arriving at a 
physically meaningful definition of particle number densities, we derive master time evolution 
equations for statistical distribution functions, which are valid to all orders in perturbation 
theory and all orders in a gradient expansion. For a scalar model, we make a loopwise truncation 
of these evolution equations, whilst still capturing fast transient behaviour, which is found to be 
dominated by energy- violating processes, leading to non-Markovian evolution of memory effects. 



1. Introduction 

The description of out-of-equilibrium many-body field-theoretic systems is of increasing relevance 
in theoretical and experimental physics at the density frontier. Examples range from the 
early Universe to the deconfined phase of QCD, the quark-gluon plasma, relevant at heavy-ion 
colliders, such as RHIC and the LHC; as well as the internal dynamics of compact astro-physical 
phenomena, such as neutron stars, and condensed matter systems. 

In [1] , the present authors introduce a new perturbative approach to non-equilibrium thermal 
quantum field theory and an alternative framework in which to derive master time evolution 
equations for macroscopic observables. In contrast to existing semi-classical approaches based 
upon the Boltzmann equation, this new approach allows the systematic incorporation of finite- 
width and off-shell effects, without the need for effective resummations. Furthermore, having 
a well-defined underlying perturbation theory that is free of pinch singularities, these time 
evolution equations may be truncated in a loopwise sense whilst retaining all orders of the time 
behaviour. Existing frameworks, based upon systems of Kadanoff-Baym equations [2j, whilst 
retaining all orders in perturbation theory, often rely upon the truncation of a gradient expansion 
in time derivatives in order to obtain calculable expressions. In this case, one necessarily makes 
assumptions as to the separation of various time-scales in these systems. In addition, one must 
generally assume a quasi-particle ansatz for the form of the propagators appearing in these 
gradient expansions. On the other hand, the loopwise-truncated evolution equations of this new 
perturbative formalism are built from non-homogeneous free propagators and time-dependent 
vertices, which together encode spatial and temporal inhomogeneity from tree-level without any 
of the aforementioned approximations. 



2. Canonical quantization 

We begin by highlighting the details of the canonical quantization of a scalar field pertinent to 
a perturbative treatment of non-equilibrium thermal field theory. 

The time-independent Schrodinger-picture field operator, denoted by a subscript S, may be 
written in the familiar plane-wave decomposition 

^s(x;t.) = /(^^ («s(p;tOe^''-'' + 4(p;t,)e-P-) , (1) 

where E{p) = y^p^ + and ag(p;tj) and as(p;ti) are the usual single-particle creation and 
annihilation operators. It is essential to emphasize that we define the Schrodinger-, Heisenberg- 
and Interaction (Dirac)-pictures to be coincident at the finite microscopic boundary time tj, i.e. 

$s(x;ti) = <Pn{ii,x;ii) = $i(ti,x;ti). (2) 

It is at this picture-independent boundary time ti that initial conditions must be specified. 
Implicit dependence on ti is marked by separation from explicit arguments with a semi-colon. 

The time- dependent interaction-picture operator is obtained via the unitary 

transformation ^i{x;ii) = e*-'^s(^'o ~ *»)<I)g(x; ?j)e~*^s(^o - ti) ^ where is the free-part of the 
Hamiltonian in the Schrodinger picture. This yields 

<^i{x;U) = (ai(p,0;t,)e-^^(P)-»e*P- + a|(p, 0; tOe^^'^P^^^e-'P- ) . (3) 

Notice that in ^ the time- dependent interaction-picture operators a\{p,i;ii) and ai(p,t;tj) are 
evaluated at the microscopic time t = 0. These operators satisfy the commutation relation 

[ai{p,i;U), a\{p',i';U)] = {27rf2E{p)6^'\p - p')e-*^(P)(* " *') , (4) 

with all other commutators vanishing, where we obtain an overall phase e~*^*-P^*-* ~ for i ^ i' . 

In quantum statistical mechanics, we are interested in the calculation of Ensemble 
Expectation Values (EEVs) of operators at a fixed microscopic time of observation tf. Such 
EEVs are obtained by taking the trace with the density operator p{tf;ti), i.e. 

(.), = Z-\t)Trp{if;U)., (5) 

where 2{t) = Tr p{tf;ti) is the partition function. Here, we have introduced the macroscopic 
time t = tf — ti, which is simply the interval of microscopic time between the specification of 
the boundary conditions and the subsequent observation of the system. 

Consider the following observable, which is the EEV of a two-point product of field operators: 

0(x,y,t/;ti) = Z-\t)TTp{if;U)^{if,^;ii)<^{if,y;U) . (6) 

It has not been necessary to specify the picture in which the operators of the RHS of ^ are to be 
interpreted. This is because these operators are evaluated at equal times. Potential observables 
built from operators evaluated at different times are picture- dependent and therefore unphysical. 
In addition, the observable O should be invariant under time translation, depending only on the 
macroscopic time t, i.e. 0{x,y,tf;ti) = O{x,y,tf — ti]0) = 0{x,y,t). Notice also that there 
are 7 independent coordinates: the spatial coordinates x and y and the macroscopic time t. It 
will later be convenient to work in terms of the central spatial coordinate X = (x + y ) /2 and 
the three-momentum q, conjugate to the relative spatial coordinate R = x — y. 
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Figure 1: The closed-time path, C = C+ U C_. The relationship between microscopic and 
macroscopic times is indicated by a dashed black arrow. 



The density operator of a time-dependent and spatially inhomogeneous background will in 
general be an intractable incoherent sum of all possible n to m multi-particle correlations, non- 
diagonal in the Fock space. We may account for our ignorance of this density operator by 
appealing to the remaining freedom in the commutation relation in ([4]). In particular, we define 



{ai{p,if;i,)al{p',if;k))t = {27rf2E{p)6'^^\p-p') + 2E'/\p)E'/\p')f{p,p',t), (7a) 
{a\{p',if;k)ai{p,if;ii))^ = 2E'^\p)E'/\p')f{p,p',t), (7b) 

where /(p, p',t) = f* (p' ,P,t) ■ The statistical distribution function f{p,p',t) is related to the 
particle number density n(q, X, t) via the Wigner transform 



n(q,X,t) = y^e^Q-^/(q + Q/2,q-Q/2,f). 



(8) 



Notice that spatial homogeneity is broken by the explicit dependence of /(p, p',i) on the two 
three- momenta p and p'. In the thermodynamic equilibrium limit, we have the correspondence: 
/(p,p',t) ^ /eq(p,p') = (27r)35{3)(p-p')/B(i?(p)) , where /b(x) = (e/^- - 1)-Ms the 
Bose-Einstein distribution function and /3 is the inverse thermodynamic temperature. 



3. Schwinger Keldysh CTP formalism 

We require a path-integral approach to generating EEVs of products of field operators. Such an 
approach is provided by the Schwinger-Keldysh CTP formalism [3l Hj . 

In order to obtain the generating functional of EEVs, we insert unitary evolution operators 
to the left and the right of the density operator in the partition function Z, yielding 



2^[p,J±,t] = Tr 



Pu{tf;ti 



Te 



«/f2 d^'x J+(a;)<I>H(a;0 



(9) 



where the spacetime hypervolume — [—t/2, t/2] X is temporally bounded. 

We may interpret this evolution as defining a closed contour C = C+ U C_ in the complex-time 
plane (t-plane, t G C), as shown in figure [l| which is the union of two anti-parallel branches: C+, 
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running from ti to tj — ie/2; and C_, running from tj — ie/2 back to ii — ie. A small imaginary 
part e = 0+ is added to separate the two, essentially coincident, branches. We may introduce 
a parametrization of this contour z{u), where u increases monotonically along C, which allows 
the definition of a path ordering operator T^. Notice that, in contrast to other interpretations 
of the CTP formalism, this contour evolves in time, with each branch having length t. 

Following the notation of [5l|6], we denote by ^±{x) = ^{x^ G C±,x) fields confined to the 
positive and negative branches of the CTP contour. We then define the doublets 



^+{x) , -<I>_(x) 



(10) 



where the CTP indices a, 5 = 1, 2 and rjab = diag(l, —1) is an SO(l, 1) 'metric' 

Inserting into ^ complete sets of eigenstates of the Heisenberg field operator, we derive 
a path-integral representation of the CTP generating functional, which depends on the path- 



ordered CTP propagator iA"^(x, y, t/; ij), written as the 2x2 matrix 



iA^\x,y,if;ii) = {Tc[<^''{x;iiW{y;ky 



AF{x,y,ij;ii) A<(x, y, t/; ti) 



(11) 



For x^, E C-^., the path-ordering Tq is equivalent to the standard time-ordering T and 
we obtain the time-ordered Feynman propagator iAp{x,y,tf;ti). On the other hand, for 



y^ £ C_, Tc is equivalent to anti-time-ordering T and we obtain the anti-time-ordered Dyson 
propagator iAY){x, y,if;ii). For x^ G C+ and y^ ^ C^, is always 'earlier' than y^, yielding the 
absolutely-ordered negative-frequency Wightman propagator iA^{x,y,if;ii). Conversely, for 
y^ E C+ and x^ S C_, we obtain the positive-frequency Wightman propagator iAy{x,y,if;ii). 

From a Legendre transform of the CTP generating functional, we derive the Cornwall-Jackiw- 
Tomboulis effective action [X\. We may then obtain the CTP Schwinger-Dyson equation: 



^abi^^V^'^f'^'>^i) = A°J, ^{x,y) + Uab{x,y,tf;ti) 



(12) 



where A^j^(x, y, t/; U) and A||^ ^(x, y) are the resummed and free inverse CTP propagators and 
^ab{x , y , t f] ti) is the CTP self-energy, analogous to (11). 



4. Non-homogeneous diagrammatics 

We consider a simple scalar theory, which comprises one heavy real scalar field <^ and one light 
pair of complex scalar fields {x^ , x), described by the Lagrangian 



X 



w?x^X 



g^x^x 



(13) 



where M = 1 GeV and m = 0.01 GeV. The ellipsis contains omitted self-interactions and the 
spacetime dependence of the fields has been suppressed. We prepare two isolated but coincident 
subsystems and S^^, both separately in thermodynamic equilibrium at the same temperature 
T = 10 GeV with the interactions switched off. The subsystem contains only the field $ and 
^y., only X- At t = 0, we turn on the interactions and the system J/' = J/'^UJy^^ re-thermalizes. 
The subsystem is taken to be infinite so that it is unperturbed by interactions with J/'^. 

The one- loop non-local $ self-energy is shown in figure [2j In particular, we draw attention 
to two features of the modified Feynman rules. Firstly, with the vertices, we associate a term 



t 

2^^ 



-smc 



(14) 



where the pi = q, ki, k2 are the four-momenta flowing into the vertex. The phase e^'^'^f , 
where qo is the energy flow external to the loop, results from the proper consideration of the 
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Figure 2: The non-local one-loop $ self-energy, ill^ ^^,((7, g', tj; tj). 
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Figure 3: The four processes contributing to the one- loop time-dependent ^ width 



Wick contraction and field-particle duality relations. Due to the finite upper and lower bounds 
on time integrals appearing in the CTP generating functional, the familiar energy-conserving 
Dirac delta function has been replaced by a sine function in (14). This violation of energy 



conservation is shown diagrammatically by the dotted line terminated in a cross and results 
from the uncertainty principle, since the observation of the system is made over a finite time 
interval. Furthermore, by virtue of this energy-violation, the perturbation series remains free of 
the pinch singularities that would otherwise result from products of delta functions with identical 
arguments at early times. Secondly, the double lines occurring in the CTP propagators of the 
loop reflect the violation of three-momentum due to the dependence on the inhomogeneous 
statistical distribution function /(p,p',t). The full set of non-homogeneous free propagators 
are listed in table [T] Together, these modified Feynman rules encode the time-dependence and 
spatial inhomogeneity of the system from tree-level. All four-momenta internal to the loop are 
integrated over and that the usual combinatorial factors apply. The CTP indices a, 6 = 1, 2 
indicate the location of the vertex on either the positive or negative branches of the CTP contour. 
For the system J^, the one-loop time-dependent <1> width is given by the following integral: 



647r3M ^ 

Ql, 02 



[ d^k smc[{qo-aiEi-a2E2)t]{l+fBiaiEi)+fBia2E2)) , (15) 

J -C/l-C/2 



±1, El 



E^Ck) 



\/k^ + and E2 



i?^(q — k). The violation of 



where ai, 02 

energy conservation, due to the sine function in (15), leads to otherwise- forbidden contributions 
from for ai, 02 = —1 (total annihilation) and ai = —02 (Landau damping). In addition, 
the kinematically-allowed phase space for normal 1 — )• 2 decays is expanded. These evanescent 
processes are shown in figure [3j For t — )• 00, we recover the known equilibrium result, since 

lim — sine r (go — ai-Ei — 02-^2) = 6(qo — aiEi — a2E2) . (16) 

t ^ 00 TT 
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Propagator Double-Momentum Representation 



Fey™ iAV,(P,P'.*7;f.) = ^^n^T^P-)'*' ''(f " P') 

+27r|2po|^/^'^(p^ - M2)/(p,y,i)e*(P«-Po)*/27r|2p^|V25(p'2 _ ^j2) 

+(-)ve-freq. iAl^^^{p,p' ,if;U) = 2^e{+{-)Po)Kp' - M^){2it)H^^){p - p') 
^^^^^"^^"^ +27r|2po|'/'<5(p2 - M2)/(p,p',t)e^(fo-Po)*/2^|2p^|i/25(p'2 _ ^2) 

Retarded i 

(Advanced) ^^r(A)(P'P') = (^^ + (_)^g)2 _ p2 _ ^2 ^^""^ ~ 
Pauli-Jordan iAO(p,p') = 2^e(po)5(p2 - M2)(27r)45(4)(p - p') 

Hadamard i^\Mrtf;U) = In 6 (p^ - M^) (27:)'' 6('\p - p') 

+27r|2po|^/^5(p^ - M2)2/(p,y,t)e^(Po-fo)t/27r|2p'o|i/25(p'2 - M^) 

Principal-part zA^(p,p') = V , \ J 2n)H^'^\p - p') 

p'^ — Al^ 

Table 1: The full complement of non-homogeneous free propagators for the scalar field <1>, where 

f{p,p',t) = e{po)e{p',)fip,p',t) + 0(-po)e(-p^)r(-p,-p',t). 

In figure 4, we plot the ratio f^^''(|q|,f) = r|jj''(|q|,t)/r^^'*(|q|,t — ^ oo) of the time-dependent 
width to its late-time equilibrium value as a function of Mt for q"^ = M^. In addition, we show 
the separate contributions of the four processes in figure [3j We note that the oscillations in 
the width have time-dependent frequencies. This non-Markovian behaviour is inherent to truly 
out-of-equilibrium systems exhibiting so-called memory effects. 

5. Master time evolution equations for particle number densities 

In order to count both on-shell and off-shell contributions systematically, we 'measure' the 
number of charges, rather than quanta of energy. This avoids any need to identify 'single- 
particle' energies by a quasi-particle approximation. We begin by relating the Noether charge 

Q{xo;ii) = ijd^x (^<^>l^{x]ii)TTl^{x]ii) - nnix;U)<^>uix]ii)^ (17) 

to a charge density operator Q(q, X,Xo;ti) via 

f d'q 



Q{Xo;ti) = j d^X Q(q,X,Xo;ti) . (18) 

By taking the equal-time EEV of Q(q, X,Xo;tj) and extracting the positive and negative 
frequency particle components, we introduce the following definition of the particle number 
density in terms of off-shell Green's functions: 



x7^t "7 2tt J (27r)4 



n(q,X,t) = Jim,2/^/^e-*«-^0(go)goiA<(Q+^,g-^,t;O), (19) 



where we have used the translational invariance of the CTP contour. 
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Figure 4: Left: the total ratio f versus Mt, for on-shell decays with |q| = 1 GeV (sohd black), 

10 GeV (blue dotted) and 100 GeV (red dashed). Right: the separate contributions to f for 
|q| = 10 GeV. The two Landau-damping contributions are equal up to numerical errors. 



By partially inverting the CTP Schwinger-Dyson equation in (12), we derive the following 
master time evolution equation for the statistical distribution function /(q + y , q — ^,t): 



5J(q+§,q-§,t) 



2 / / dft, dQo q . Q (g + I, g - f , t; 0) 



2tt 2tt 



(20) 



2 ) y 2 ' 



t;0) 



ffdqodQo iQ^t 



J J 2tt 2tt 
where we have introduced 



0{qo 



+ 



2 ' ' 



Q 
2 ' 



t;0) + ^* 



2 ) y 2 ' 



t;0) 



^{q+i,q-^,t;0) = - iUp{q + Q,k,t;0) iA<{k,q - Q,t;0) , (21a) 



+ 



2 '' 



',t;0) 



1 r d^k 

27 (2^ 



m>(g+ f ,/c,t;0) iA<{k,q- f ,t;0) 



in<(g+ 2,A;,t;0) ( iA>(A:,g- f ,t;0) - 2iAp{k,q 



2 

|,t;0) 



(21b) 



It is important to stress here that (20) provides a self-consistent time evolution equation for / 



valid to all orders in perturbation theory and to all orders in gradient expansion. The terms on 



the LHS of (20) may be associated with the total derivative in the phase-space (X, p), which 
appears in the classical Boltzmann transport equation. The ^ terms on the LHS of (20) are 



the force terms, generated by the potential due to the dispersive part of the self-energy and the 



^ terms n the RHS of (20) are the collision terms. 
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Truncating (20) to leading order in a loopwise sense, we obtain for our simple scalar theory, 
the following time evolution equation for the <I> statistical distribution function: 



- E 

2 ^ 

t . 

X — smc 
2tt 



d^k 



(27r)3 2£;$(q) 2^^(k) 2£;x(q - k) 
[(a^$(q) - ai^x(k) - a2^x(q " k))V2 



|7r + 2Si (^aSci,(q) + ai^x(k) +a2-Ex(q-k))t/2 | 
{[^(-a) + /$(|qM)] [e(ai)(l + /x(|k|,t))+^(-ai)4^(|k|,t) 

X [0(a2)(l + /^(|q - k|,t)) + ^(-a2)/x(|q - k|,t)] 
[^(a) + /*(|q|,t)] [^(ai)/x(|k|,t) + 0(-ai)(l + /^(|k|,t))] 

X [0(a2)/^(|q - k|, i) + ^(-a2)(l + /x(|q " k|, t))] }, 



(22) 



where a, ai, a2 = ±1- The second and third lines of (22) encode the early-time violation 



of energy conservation. This leads to the non-Markovian evolution of memory effects and 
evanescent contributions from otherwise kinematically-disallowed processes. Replacing these 
lines by an energy-conserving delta function, we recover the semi-classical Boltzmann transport 
equation. However, given the equilibrium initial conditions of our model, the artificial imposition 
of energy conservation along with the properties of the Bose-Einstein distribution ensure that 



the RHS of (22) is zero for all times. Thus, the semi-classical Boltzmann equation cannot 



describe the re-thermalization of our simple model. This is true also for gradient expansions of 
Kadanoff-Baym equations when truncated to zeroth order in time-derivatives. Hence, it is only 
when energy-violating effects are systematically included, as in this new perturbative approach, 
that the dynamics of this re-thermalization are captured. 



6. Conclusions 

We have obtained master time evolution equations for particle number densities that are valid to 
all orders in perturbation theory and to all orders in gradient expansion. The loopwise truncation 
of these time evolution equations remains valid to all orders in gradient expansion, capturing 
the evolution on all time scales, including the transient dynamics. This prompt behaviour 
is dominated by energy-violating processes that lead to non-Markovian evolution of memory 
effects. The underlying perturbation series are built from non- homogeneous free propagators 
and explicitly time-dependent vertices. Due to the systematic treatment of finite boundary and 
observation times, these diagrammatic series remain free of pinch singularities. 



Acknowledgments 

The work of PM and AP is supported in part by the Lancaster-Manchester-Sheffield Consortium 
for Fundamental Physics under STFC grant ST/J000418/1. AP also acknowledges partial 
support by an IPPP associateship from Durham University. 



References 

[1] Millington P and Pilaftsis A 2012 {Preprint 1211 .3152 ) 

[2] Kadanoff L P and Baym G 1989 Quantum statistical mechanics (Addison- Wesley) 

[3] Schwinger J S 1961 J.Math.Phys. 2 407-432 

[4] Keldysh L 1964 Zh.Eksp.Teor.Fiz. 47 1515-1527 

[5] Calzetta E and Hu B L 1987 Phys.Rev. D35 495-509 

[6] Calzetta E and Hu B L 1988 Phys.Rev. D37 2878-2900 

[7] Cornwall J M, Jackiw R and Tomboulis E 1974 Phys.Rev. DIO 2428-2445 



8 



